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A CLASS OF NONLINEAR OPTIMIZATION PROBLEMS WITH 

POTENTIALS 

JIAKUN LIU 



o 

^SJ ' Abstract. In this paper, we introduce a class of nonlinear optimization problems. Under 

mild assumptions, we obtain the existence of potential functions and show that the poten- 
cy ■ tial function is a generalized solution of a Monge- Ampere type equation. We also present 
some interesting applications in optimal transportation and geometric optics problems. 



1. Introduction 

In a previous paper [16], we proved that the reflector problem in the near field case is a 
nonlinear optimization problem. In this paper we introduce a class of nonlinear optimization 
problems with potentials and prove the existence of optimal mappings. These problems 
are natural extensions of the optimal transportation problem, which is a class of linear 
optimization problems extensively studied in recent years. In particular, we give more 
examples of nonlinear optimization problems arising in reflector and refractor problems. 

•/^ ' Our results, together with jl2| 125]. show that the light reflection and refraction processes 

^v^ , can be formulated to nonlinear optimization problems. 

cn . Let's consider a class of optimization problems of maximizing the following functional 

^ ' (1-1) Iiu,v) := / F{x,y,u{x),v{y))d'y, 

Juxv 

where U, V are two bounded domains in M", F is a function in ]R^"+^, dj is a measure on 

K> ■ U X V, which has dx, dy as its marginals, and (u, v) is a pair of functions contained in the 

b ■ 

C^ ' constraint set 

(1.2) K = {{u,v)e C{U) X C{V) : 0(x, y,u,v) < in U x V} , 

where (j) is the nonlinear constraint function. More generally, we may also consider U, V 
subsets of a Riemannian manifold {M , g) . Some interesting examples and applications of 
the problem p.ip - p.2p are included in Section [5l 

Write F = F{x,y,t,s) and cj) = (j){x,y,t,s), where x,y,t,s are independent variables. 
Use the subscripts to denote the partial derivatives, i.e., Fx^ = dF/dxi, (pt = d(j)/dt, etc. 
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We always assume the following hypotheses: The function F is C^ smooth and monotone 
increasing in t,s, namely 

(1.3) Ft>0, Fs>0 in U xV xR^, 

and the constraint cp = (j){x, y, t, s) is C^ smooth and strictly increasing in t, s. Note that 
by the implicit function theorem, the constraint <j){x, y, u,v) = can always be written as 

(t>(.x,y,-ip{x,y,v),v) =0, 

for some function cp = ip{x,y,s), which is C^ smooth in s. Since </> is strictly increasing in 
t, without loss of generality, the constraint <j) < can be written as 

(1.4) (l){x,y,u,v) =u + ip{x,y,v) <0. 

Correspondingly, the function ip = (p{x, y, s) is strictly increasing in s. We assume further 
that there exists a constant ^o > such that 

(1.5) V's > 6*0 in U xV xR. 

In order to state our main results, we introduce the following conditions and define some 
terminology. 

(HI) For each xq € U, for any (p, t) G M" x M, there is at most one pair {y, s) G M" x M 
such that 

{ipx,ip){xo,y,s) = -{p,t), 

namely, ipxixo,y,s) = -p and ip{xo,y,s) = -t. 
(H2) For any {x,y,s) eU xV x M, 

det ( ipxy (fxs (Sxpy) / 0. 

Definition 1.1. A pair {u,v) £ K is called a dual pair with respect to (j) if 

(1.6) u{x) = sup{i : (l){x,y,t,v{y)) < 0, Vy G V}, 
v{y) = sup{s : (/>(x, y, u{x), s) < 0, Vx G U}. 

If furthermore I{u, v) = sup;^ /, {u, v) is called a dual maximizing pair of I, and the 
functions u, v are called potential functions in the nonlinear optimization (jl.ip - (ll.2p . 

Definition 1.2. Let {u,v) £ K be a dual maximizing pair of F A mapping T : U ^ V is 
called an optimal mapping associated to (u, v) if T solves the equation 

(1.7) (/)(x,Ta:;,n(x),'y(rx)) = 
almost everywhere on U . 
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The following theorem gives the existence of potential functions and optimal mappings 
in the nonlinear optimization problem (|l.ip - (jl.2p . Their regularity properties will be inves- 
tigated in a subsequent paper [T7] . 

Theorem 1.1. In addition to hypotheses (jl.3p - ()1.5p and (H1)-(H2), we assume the fol- 
lowing mixed condition on F and 93; for any (u, v) E K , there holds 

(1.8) / \-Ft{x,y,u{x),v{y)) + —{x,y,u{x),v{y))\d-i = Q. 
Juxv I V^s J 

Then there exists a dual maximizing pair (n, v) £ K of I, and a unique optimal mapping T 
associated to {u, v) uniquely determined almost everywhere on U . 

Let us consider a special case, which has many practical applications, that 

(1.9) F{x, y, u, v) = u{x)f{x) + Lp{x, y, v)g{y), 

where / E L^{U), g E L^{V), and if is the constraint in (|1.4p . By the monotonicity 
(jl.3p - (jl.5p . we assume that f,g are positive functions. The condition (jl.Sp implies 



0=/ {-f{x)+g{y)}dj. 
JUxV 

Since d^ has dx, dy as its marginals, it is equivalent to the balance condition 

(1.10) [ f{x)dx = I g{y)dy. 

Ju Jv 

That is to say, in the case of ()1.9p the condition (jl.Sp , or equivalently (jl.lOp , is independent 
of the functions u^v. 

As a consequence of Theorem ll.H we have 



Corollary 1.1. Assume that the function F is given by (jl.9p for some positive functions 
f E L^{U),g E L^{V) satisfying the balance condition (jl.lOp . Assume that the constraint if 
satisfies ()1.5p and (H1)-(H2). Then the same conclusions in Theorem \l.l\ hold. 

We point out that for the nonlinear optimization ()l.ip - (jl.2p . there is generally no unique- 
ness for maximizing pair {u,v), see Remark 12.11 Under the following convexity hypothesis 

(He) if is convex in s, and F is strictly concave in (t, s)-variable, 

one can easily obtain the uniqueness of maximizing pair. However, this hypothesis (He) 
never holds in the case of (jl.9p . The proof of uniqueness under (He) is straight-forward. 
Suppose {u, v) £ K and {u, v) £ K are two different maximizing pair, namely 

(1.11) I{u,v) = I{u,v) = sup I =:!*. 

K 

Set u = 2(u + u) and v = ^{v + v). Since ^ is convex in s, we have 

(1.12) 0(x, y, u,v) <- ((^(x, y, u, v) + (p{x, y, u, v)) < 0, \/{x, y) £ U x V, 
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which imphes that (n, v) € K. On the other hand, since F is strictly concave in (t, s), 

(1.13) I{u,v)>^{I{u,v) + I{u,v)) = I*, 
which is a contradiction to (jl.lip . 

It is weh known that many constrained nonlinear optimization problems can be solved by 
Lagrangian dual methods (see, for example, [lO]), where the convexity plays a crucial role. 
For the nonlinear optimization (|l.ip - (jl.2p under conditions (|1.18p - (|1.19p . which is weaker 
than the hypothesis (He), we can show that there is no duality gap and there exists at least 
one Lagrange multiplier. This enables us to use the Lagrangian duality theory to study the 
maximization of the functional /, see Section HI 

Now, let us introduce some terminology in Lagrangian duality. More details are con- 
tained in Section m Denote X := C{U) x CiV). The nonlinear optimization (|l.ip - (|1.2p is 
equivalent to the following primal problem: 

maximize I{u,v), 

(1.14) subject to {u,v) G X, ip{u,v) := inf —(p{x,y,u,v) > 0. 

UxV 

We always assume that (j) has the form (jl.4p . An element (u,v) G X is called feasible if 
ip{u, v) > 0, or equivalently, {u, v) G K. 

Define the Lagrangian function L : X x M — ;• M to be 

(1.15) L{u,v, fi) = I{u,v) + fj.ip{u,v), 
where /z G M. The dual functional J is defined by 

(1.16) J{fi) = sup L{u,v,fi), 

(u,v)ex 

and the dual problem is given by 

minimize J(/i) 

(1.17) 

subject to /z > 0. 

Regardless of the functional / and the constraint (p of the primal problem, the dual problem 
has a very nice convexity property, as shown in Lemma 14. 2i In the language of nonlinear 
programming [3l[T0], when inf^>o J(/i) = sup(„ ,^)g^ /(ti,f), we say that there is no duality 
gap, otherwise, there is duality gap. 

Theorem 1.2. Assume that the function F in (jl.ip is concave in {t,s), namely for any 
{x,y,t,s) eU xV xR"^, 

(1.18) ffe,.,„F:=(|;5 ^p ) < 0. 

and the constraint if in (jl.4p is convex in s, namely for any {x,y,s) (i U x V x M, 

(1.19) dss^ > 0. 
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Suppose that there exists a pair (u, v) £ X such that 

(1.20) ^(m,u)>0. 

Then there is no duality gap between the primal problem (J1.14p and dual problem (J1.17p . 
and there exists at least one Lagrange multiplier, (see Definition\4.1\). 



Note that in the special case (jl.Op . the convexity assumptions (jl.lSp and (|1.19p imply 
that (p = (p{x, y, v) is linear in v, namely 

ip{x, y, v) = Co{x, y) + ci{x, y)v, 

for some functions cq,ci, where ci > 6q > in U x V hy the monotonicity (ll.Sp . When 
ci = 1, it is an optimal transportation problem with the associated cost function — cq, see 
Example 15.11 

This paper is organized as follows: In Section [2] we show various properties of potential 
functions and optimal mappings under appropriate assumptions. In particular, we show the 
optimal mapping T is measure preserving in the sense of ()2.7p and the potential function 
satisfies a Monge- Ampere type equation. Theorem II. II is proved. In Section [3] we introduce 
a notion of generalized solutions and show that a potential function is a generalized solu- 
tion, from which the existence of generalized solutions follows, see Theorem 13.11 Section 
m contains the Lagrangian duality theory. It provides a useful tool to obtain the existence 
of maximizer. Theorem 11.21 is proved. In Section [5] we present some examples and ap- 
plications of nonlinear optimization (ll.ip - (ll.2p . In particular, we derive the equations in 
geometric optics problems from the corresponding constraints, instead of using the reflection 
or refraction law. 

Acknowledgements. The author would like to thank Xu-Jia Wang for proposing the prob- 
lem to him. He would also like to thank Neil Trudinger for discussions on the topic at 
the conference at Hefei, China in August 2011. This work was supported by the Simons 
Foundation. 



2. Existence of maximizers 

In this section, we study the existence of potential functions and optimal mappings in 
the nonlinear optimization problem (jl.ip ~^ (jl.2p . 



Lemma 2.1. Assume that the function F and the constraint (/) satisfy the hypotheses (|1.3p - 
(jl.Sp and (jl.Sp . Then the objective functional I in (jl.ip has a maximizing pair {u,v) £ K, 
where K is the constraint set (jl.2p . 
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Proof. Given any pair {u, v) G K, we claim that I{u, v) does not decrease if v is replaced 
by 

(2.1) v*{y) = sup{s : <j){x , y , u{x) , s) < 0, Vx G U}. 

In fact, by the continuity of (p and u, for each y G V there is some x G U such that 

4'{x,y,u{x),v*{y)) = > cl){x,y,u{x),v(y)), 

since (n, -y) G K. By ()1.4p - ()1.5p . v* > v. Furthermore, (j){x , y , u{x) , v* (y)) < for all 
(x,y) £U xV, so {u,v*) £ K. 

Since f * > f, by (jl.Sp we have 



/(n,u*) >Iiu,v). 
Similarly, if we define 

(2.2) u*{x) = sup{t : (P{x,y,t,v*{y)) < 0, Vy G F}, 

then {u*,v*) G K and 

/(n*,?;*) >I{u,v*) >I{u,v). 

Thus we do not decrease I{u,v) by replacing {u,v) by (n*,t'*). The claim is proved. 

Define Kc^ = K r]{u > Co}, where Cq is a constant, which may be chosen negative and 
sufficiently small in the following context. We show that u* and v* are uniformly bounded 
if {u,v) G Kco- Since t;* > i;,n > Cq, by (jl.4p - (jl.5p we have for each y gV, s := v*{y), 

Cq + if{x, y, s) < u{x) + if{x, y, s) < 0, for all x £ U. 

Then by (jl.Sp again, there exists a constant Ci, such that s < Ci. This implies that 

(2.3) v<v* <Ci, 

we may choose Ci such that supy v* = Ci. By a similar argument, there is another constant 
Co depending on ip and Ci such that inf[/ u* = Cq. The constant Cq > Cq, since u* > u in 
U, and so (n*,f*) G Kcq- 

We next deduce the lower bound of v* and the upper bound of u* by showing that u* and 
V* are locally Lipschitz. Consider two points in U, xi ^ X2 and \xi — X2\ < £ sufficiently 
small. There are two points yi,y2 ^V such that 

cl){xi,yi,u*{xi),v*{yi)) = 0, 
(I>{x2,y2,u*{x2),v*{y2)) = 0. 



NONLINEAR OPTIMIZATION WITH POTENTIALS 7 

Then by (jl.4p , we have 

=(p{x2 ,y2,u* (x2) , V* (y2)) - H^i ,y2,u* (xi ),V*{y2)) 

+ 4>{xi,y2,u* (xi),v* {y2)) - (l){xi,yi,u*{xi),v*{yi)) 

=u*{x2) -u*{xi) -ipx{x,y2,v*{y2)) ■ {x2 - Xi) 

+ 4>{xi,y2,u*(xi),v*{y2)), 

where x = 6xi + (1 — 0)x2 for some 6 £ (0, 1). Noting that 4>{xi,y2,u*{xi),v*{y2)) < 0, we 
have 

U*{X2) - U*{xi) > -C2\X2 - Xl\, 

where the constant C2 = sup{\(px\ + \^y\)- 

On the other hand, replacmg 4>{xi,y2,u*{xi),v*{y2)) by (I){x2,yi,u*{x2),v*{yi)) in the 
above calculation, we have 

U*{X2) - U*{xi) < C2IX2 - Xi|. 

Therefore, the Lipschitz constant of u* on U is controlled by 

(2.4) \\u*\\Lip^u)<C2- 

By switching x and y in the above argument, we can obtain the Lipschitz continuity of v* 
on V, 

\v (2/2) - V [yi)\ < '' \y2 - yi\ 

mt(/?s 

< C20Q^\y2 -yi\, 

where ^0 is the constant in (|1.5|) . 2/1,2/2 are two distinct points in V. This inequality 
implies that ||f*||2,ip(y) < (^2^^ ■ Hence, we have u* < Cq + C2diam(C/) and v* > Ci — 
C2eQ^d[am{V). 

We conclude, therefore, that any pair [u, v) G Kc^ may be replaced by a bounded, 
Lipschitz pair {u*,v*) G Kqq without decreasing /. We now choose a maximizing sequence 
{iuk,Vk)} C Kco such that 

I{uk,Vk)^ sup I{u,v). 

By the above considerations we may assume that each {uk,V}S) is a bounded, uniformly 
Lipschitz pair, uniformly with respect to fc, so there is a subsequence converging uniformly 
to a bounded, Lipschitz, maximizing pair {u^v) E -f^Co- 
Last, we show that when Cq < sufficiently small, 

sup I{u,v) = sup I{u,v), 
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or equivalently, sup;^^ I is independent of Cq. By definition, one has sup^^ _^I > 
sup^^ I. So it suffices to sliow tfie reverse inequality. Let {u^v) G i^Co-i be a maxi- 
mizer such that /(n, v) = supj^ _ /, and {xk}k=i,--- ,N be a set of points in U . For a smaU 
constant e > 0, define 

u inU -UNBeixk), 



n + 2 in UAri?e(xfc) 

Note that we may replace u by its mollification Uh = Ph * u, where ph is the standard 
mollifier function [U]. For simplicity, we assume u continuous in the sense that for h > 
sufficiently small, 

I{uh,v) = I{u,v) + 0{Ne''). 
Define 

^*{y) = sup{s : (p{x,y,u{x),s) < 0, yx£U}, 
u*{x) = sup{t : (l){x,y,t,v*{y)) <0, \/y€V}. 

Since the constraint function 99 is C^ smooth in s and by (jl.4p - (jl.5p . except a set E C U 
and a set E' C V of measure |ii;| = lii^'l = 0{Ne'^), 

V* = v- — + 0{6) mV\E', 

u* = U + 2 + 0{5) mU\E, 

where 5 := minj^j{dist(xj,Xj)}. Therefore, by (jl.Sp and the mean value theorem we have 

Iiu*,v*) = Iiu,v) + 2 [ I Ft - —] dj + 0(6) + 0{Ne'') 

JiU\E)x{V\E') I V^s ) 

>I{u,v)-C6-CNe''. 

As {u, v) £ Kcq-i, we may assume that inf[/ u = Cq — I. Otherwise, one has inf^/ u = Cq — tq 
for some constant tq < 1. This implies that sup^^ _^I = sup^^ I, namely sup^^ / is 
independent of Co, and the proof is finished. By the definition, 6 will become small if the 
number of points A^ is sufficiently large so that we have (u* ,v*) £ Kqq and 

sup/>/(r,r)> sup /-CJ-CiVe". 

Kcq ^Co-1 

Then, choosing e > sufficiently small, we have 

sup / < sup /, 

Kco-i Kcg 

by letting 5 — )■ 0, e — )• 0, which implies that sup^^ / is independent of Co, and the proof is 
finished. D 

Remark 2.1. From the proof of Lemma |2.H we conclude that there exist infinitely many 
maximizing pairs. In fact, if {u,v) is a maximizer and Cq = infj/u, then there is another 
maximizer in Kcq+i, which is different from {u,v). 
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Lemma 2.2. Let (u, v) £ K be a dual maximizing pair of I in Lemma \2.1[ Assum,e ip 
satisfies condition (HI). There exists an associated mapping T with respect to {u,v), which 
solves equation (jl.7p and is uniquely determined almost everywhere on U. 



Proof. Since u satisfies (I1.6P and v, ip are continuous, for each x G C/, there exists some 
y =: T(x) G V such that 

(2.5) u{x) + ip{x,y,v{y)) = 0, 

u{x') + ip{x',y,v{y)) < 0, 

for any other x' £ U. Note that u is differentiable almost everywhere since u is Lipschitz 
continuous. Let x £ U he a differentiable point of u, by differentiation we have 

(2.6) ipx{x,y,v) + Du{x) = 0. 

Therefore, for the fixed x £ U, setting t = u{x) and p = Du{x) one can see that 

(px{x, y, v) = -p, and (p{x, y, v) = -t. 

From the condition (HI), we then obtain the uniqueness of y = T{x). Since u is differentiable 
almost everywhere on U, the mapping T is determined almost everywhere on U. D 

Remark 2.2. Li the condition (HI), we assume that {ipx,^){x, •, •) is one-to-one in the whole 
space M" x M, for each x £ U. This is only for simplicity. We may allow that the constraint 
(p is defined in a proper subset V( xl, where 1/1 C M" x M" and X C M. Denote the projections 
Ux = {y£W^ : {x,y) £ U], Uy = {x £ W : {x,y) £ U], let U = [jUy, V = (J^Cx- In this 
case we replace (HI) by assuming that: for any (x,y) £ U, there exists an open interval 
I{x,y) C X, such that {ipx,(p){x,-,-) is one-to-one in y £ Ux,s £ I{x,y), for each x £ U. 
Accordingly the condition (H2) can be restricted to ZY x X, namely for any (x, y) £ U, 
s £ X(x, y), we have det {ipxy — (pj^^pxs ^ ^y) 7^ 0. 

Let (n, v) £ K he a dual maximizing pair of /, and S : U ^ V. We say the mapping 5 
is measure preserving if 

(2.7) 0=/ {-Ftips{x,Sx,v{Sx))h{Sx) + Fs{x,y,u{x),v{y))h{y)}dj, 

JuxV 

for any h £ C{V). 

Lemma 2.3. Assume that the balance condition (jl.Sp holds. Let T be the optimal mapping 
obtained in Lemma \2.2\, associated with a dual maximizing pair {u,v). Then T is measure 
preserving in the sense of ()2.7p . 



Proof. Let h £ C(y) and |e| < 1 sufficiently small. Define 
(2.8) v,{y) = v{y) + eh{y) 
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and 

(2.9) u,{x) = sup{t : t + ip{x, y, v,{y)) < 0, Vy G V]. 

Then {ue,Ve) S K and (uq^vq) = {u,v). 

Since {u,v) satisfies (|1.6p . by Lemma [2^2] for every x £ U the supremum (jl.6p is attained 
at point yo = T{x). We claim that at these points we have 

(2.10) Ue{x) - u{x) = -eipsh{T{x)) + o(e). 

To prove ([230]) . first we show that LHS < RHS. 

= u{x) +ip{x,yo,v{yo)) 
= u{x) + ip{x, yo, Ve{yo) - eh{yo)) 
= u{x) + ip{x, yo,Ve{yo)) - e0s/i(yo) + o{e) 
= u{x) - n,(x) + (/)(x, yo,Ue{x),Ve{yo)) - e^sHyo) + o(e) 
< u{x) - n,(x) - eipsh{yo) + o(e). 

Thus we have 

Ue{x) - u{x) < -eipsh{yo) + o(e). 

To show LHS > RHS we use the fact that for any such x £ U there are points y^ £ V 
such that the supremum in (|2.9p is attained. Thus 



Then we have 



> u{x) + ip{x,y^,v{ye)) 
= u{x) + ^{x,ye,Ve{ye) - ^h{y^)) 
= u{x) + ip{x, y^,Ve{ye)) - eipsh{y^) + o(e) 
= u{x) - Ue{x) + (j){x, ye, u^{x) , v^{y^)) - e(psh{ye) + o(e) 
= u{x) - Ue{x) - eLpsh{ye) + o(e). 

u^{x) - u{x) > -€Lpsh{ye) + o(e) 

= -eifsHyo) + eifs (/i(yo) - HVe)) + o{e). 

Since the supremum in ()1.6p is attained at yo, we have y^ — )■ yo as e — t- 0, and therefore, 
since h G C{V), 

eips{h{yo) -h{ye)) = o(e). 
This implies that LHS > RHS, and (l2T0]l follows. 
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Next, since {u,v) = {uo,vo) maximizes / among all pairs in K, we obtain 
= lijn I{u„v,)-I{u,v) 

^ ^.^ f F{X, y, Ue,Ve) - F{x, y, u, v) ^ 



{-Ft^sh{T{x)) + FMy)}d-f. 
UxV 



D 



In the special case of ()1.9p . by (jl.Sp and the arbitrariness of test function h G C{V) in 
()2.7p we have 

(2.11) / f{x)h{Sx)dx = [ g{y)h{y)dy, 

Ju Jv 

for any h G C{V), where f,g are positive functions in L^[U), L^{V), respectively, and satisfy 
the balance condition (|1.10p . Moreover, we have the following 

Corollary 2.1. Assume the function F is given by ()1.9p with / > 0, G L^{U), g > 0,£ 
L^{V) satisfying (jl.lOp . The optimal mapping T in Lemma \2.^ is measure preserving in 
the sense of ()2.1ip . IfT is continuous differentiable, then 

f 



(2.12) \detDT\ 



goT- 



Proof. From the proof of Lemma [2. 31 it is obvious that T is measure preserving in the sense 
of ()2.1ip . When T is C^ smooth, by the formula of change of coordinates. 



f{x)h{Tx)dx= / g{Tx)h{Tx)\detDT\dx, 
u Ju 

for any h G CiV). Hence the Jacobian of DT satisfies (j2.12p . D 

As a consequence of Lemma 12.31 and Corollary 12.11 we derive the equation satisfied by 
the potential function u as follows. At this stage, let us assume all the functions are smooth 
enough, say at least C^, so that we can do the differentiations. 

Let {u, v) £ K he the dual maximizing pair of /, and T be the associated optimal 
mapping. By Definition 11.21 and (|2.6p we have 

ipx{x, Tx, v{Tx)) + Du{x) = 0, 

in U. Differentiating with respect to x, we then get 

(2.13) = ip^^ + ipxyDT + {ipxs 'S) Dv)DT + D'^u, 

where each side is regarded as an n x n matrix valued at (x,y), y = T(x). 
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In order to eliminate Dv in (j2.13p . we note that for xq ^ U equality ()1.7p holds at 
2/0 = T{xq), and for other y' ^V 

u{xq) + ip{xo,y',v{y')) < 0, 

since {u,v) £ K. Thus, at (xo,yo) there holds 

0. 



dip ^ 

dy 



By the assumption (jl.Sp . (/J^ > we thus get 

(2.14) 






Combining (I2.13P and (I2.14p we have the equation 



(2.15) 



\D u + Lp^ 






\DT\ 



In the case of ()1.9p - (|1.10p . by Corollary 12. II we obtain the equation 



(2.16) 



|det [D^u + (^2 



det 



^xy ^xs ® (fy 



f 



goT 



Under the condition (H2) and f,g > one see that (j2.16p is a non-degenerate Monge- 
Ampere type equation [U [11] . When the optimal mapping T is a diffeomorphism from U 
to V, we have the natural boundary condition 

(2.17) T{U) = V. 

Similarly, one can also derive the dual PDE for the dual potential v. 

When ip{x, y,v) = x-y, then (|2.16p is equivalent to the standard Monge- Ampere equation 

det D u = h, 

with the boundary condition 

Du{U) = V, 

where h = f/g. When (p{x, y,v) = v — c{x, y) for some function c : [/ x y — )• M, we have 
the optimal transportation equation, see Example 15. H 

/ 



det [D'^u - Die] = \det Dlyc\ 



goT' 



with the boundary condition 

T{U) = V. 

It is well-known that the regularity of equation ()2.16p depends crucially on the structure 
of constraint function ip, as in |19( I24j . In the next section, we introduce a notion of a 
generalized solution and show that if {u, v) is a dual maximizing pair of / over K, then 
the potential u is a generalized solution of ()2.16p . The regularity property of u will be 
investigated in a separate paper [T7] . 
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3. Generalized solution 

In this section, we introduce a notion of generalized solutions of (I2.16P and show that 
the potential function is a generalized solution. Let Lp be the constraint in (|1.4p . First we 
introduce the (/3-concavity for functions, which is an extension of the c-concavity in optimal 
transportation, see [HI [19] . 

Definition 3.1. A ip-support function of u at xq is a function of the form. (p{x,yo,so), 
where yo £ I^"? '^i^'d sq £ M is a constant such that 

u{xo) + ip{xo,yo,so) = 0, 

(3.1) u{x) + ip{x,yo,so) < Vx G [/. 

A continuous function u defined on U is (p-concave if for any point xq G U , there exists 
a if-support function at xq. 

By definition, the potential function u is 99-concave with yo gV,S() = f (yo)- In the special 
case when (p{x, y,s) = s — x ■ y, the notion of c^-concavity coincides with that of concavity, 
and the graph of a (^-support function is a support hyperplane. 

Recall that ip is derived from the constraint function </>(x, y, u, v) by the strict monotonic- 
ity in u. Since <j) is also strictly increasing in v, the constraint (|1.4p can also be written 
as 

(3.2) 4>{x,y,u,v) =v + ip*{x,y,u) <0, 

for a function ip* = (p*{x, y, t) strictly increasing in t. The function ip* = (p*{x, y, t) is called 
dual constraint function of ip in the sense of 

-^{x,y,-'P*{x,y,t)) = t, 

-^*{x,y,-ip{x,y,s)) = s, 

for all {x, y) & U X V. By differentiating, we have 

(3.3 iPt = — , ip^ = — , iPy = — . 

For the dual constraint (/?*, we assume the following condition analogous to (HI): 

(HI*) For each 7/0 G y ■, for any {q, s) G M" x M, there is at most one pair {x, t) G M" x M 
such that 

(v^V?*) {x,yo,t) = -{q,s), 
namely, ipl{x,yo,t) = -q and ip*{x,yo,t) = -s. 

The c/p-concavity in Definition 13.11 and (|3.ip - (j3.3p are generalizations of c-concavity and 
c-duality in optimal transportation, where 

'P{x,y,s) = s-c{x,y), ip*{x,y,t) =t-c{x,y), 
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for a cost function c{x,y). The condition (Hl)-(Hl*) is the counterpart of the condition 
(Al) assumed on the cost function c{x,y) in [19]. Note that from ()3.2p and (j3.3p . we can 
directly derive ()2.14p for a dual pair of functions u, v. 

Similarly, by switching x and y, U and V, one can also introduce the notion of (p*- 
concavity for the function v. From Definition 11.11 and (jS.ip . when {u,v) S i^ is a dual pair, 
u is naturally 99-concave and v is c^*-concave. 

Let uhe a (/3-concave function in U. We define a set-valued mapping Tu = Tu,^ -.U^V. 
For any xq G U, let Tu{xq) denote the set of points yo such that y5(x,yo,so) is a 99-support 
function of u at xq for some constant sq. For any subset E C U, we denote Tu{E) = 

UxeETuix). 

If u is C^ smooth, by condition (HI) (yojSo) is uniquely determined by {Du{xo),u{xq)), 
and Tu is single valued. In this paper we call the mapping Tu the ip-normal mapping of u. 
Similarly we can define the 99*-normal mapping for 93*-concave functions. In particular, if 
{u,v) £ K IS a dual pair, we see that y G Tu^y,{x) if and only if x G Tu^^p*{y). 

Remark 3.1. As the constraint function ip is smooth, any 93-concave function u is semi- 
concave, namely there exists a constant C such that u{x) — C|xp is concave. It follows that 
u is twice differentiable almost everywhere and Tu{x) is a singleton for almost all x £ U. 

Lemma 3.1. Let {u,v) £ K be a dual maxim,izing pair of T Assum,e that the constraint 
ip* satisfies condition (HI*). Let 

Y = Yu = {y £ V \3 xi ^ X2 £ U such that y £ Tu{xi) n Tu{x2)} ■ 

Then Y has Lebesgue measure zero. 

Proof. If y £ Tu{xi) n Tu{x2), we have xi,X2 £ Tu^^*{y). From the proof of Lemma [27n v 
is almost everywhere differentiable. Assume y is a differentiable point, then by definition 

ip*{xi,y,u{xi)) = -v{y), 
ip*y{xi,y,u{xi)) = -Dv{y), 

for i = 1,2. If xi 7^ X2, this is a contradiction to (HI*). D 

We now define a measure /x = iiu,g in U, where g £ L^(V) is the positive measurable 
function in (fLgjl - ifrTOj) . Set 5 = in W - V. For any Borel set E CU, define 



(3.4) ^{E) = / giy)dy. 

■JTu{E) 

It follows from Lemma |3. II that fi is a Radon measure, and satisfies the following regularity 
properties: 

H{E) = mf{n{D) : E C D C U, I? open} 
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for all Borel sets E C U, and 

/u(D) = sup{/i(i^) : K C D, incompact} 

for all open sets D C U. For further discussion of the measure /i and its stability property, 
see 121 El in 19]. 

Definition 3.2. A ip-concave function u is called a generalized solution of (j2.16p if ^u,g = 
fdx in the sense of measure, that is for any Borel set E <ZU , 



(3.5) / / = / g. 

JE JTu{E) 

Note that since we extended 5 = to M — T/, the boundary condition (j2.17p is a conse- 
quence of the mass balance condition (|1.10|) in the sense of |y — Tu{U)\ = 0. 

The next result shows that a potential function u is a generalized solution. The existence 
of potential in Lemma F'2 . 1 1 then implies the existence of generalized solutions. But in general 
we do not have the uniqueness, see Remark 12.11 

Theorem 3.1. Let {u,v) £ K be a dual maximizing pair of I. Then u is a generalized 
solution of (j2.16p . 

Proof. Let {u,v) £ K he a dual maximizing pair of /. Then by ()1.6p . u is (/^-convex and v 
is Lf* -convex with respect to each other. By Lemma 12.21 the optimal mapping T associated 
to {u,v), as determined by ()2.6p . is equal to the mapping r„^^ almost everywhere on U. 
By Corollary 12.11 T is measure preserving in the sense of (|2.1ip . Hence u is a generalized 
solution of l^7rU\i . Assumption pTTUl) imphes that (pTTp holds. D 

4. Lagrangian duality 

In this section, we study the dual problem (|1.17p of the constrained nonlinear optimization 
(jl.ip - (jl.2p . and prove Theorem 11.21 Recall that the Lagrangian function L is defined in 
(J1.15p . where the constraint ^ is given in (|1.14p . Denote by /* the optimal value of the 
primal problem (J1.14p . namely 

/* = sup I{u,v). 

Definition 4.1. A factor fj,* is called a Lagrange multiplier for the primal problem if fj,* > 0, 
and 

I* = su];){L{u,v, fi*) : {u,v) G X}. 
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Lemma 4.1. Let /i* be a Lagrange multiplier. Then {u*,v*) is a global maximum of the 
primal problem if and only if (n*, v*) is feasible and 

(4.1) {u*,v*) = arg max L{u,v,n*), 

{u,v)£X 

(4.2) fi*il;{u*,v*) = 0. 

Proof. If {u*,v*) is a global maximum of the primal problem, then {u*,v*) is feasible and 
furthermore, 

r =L{u*,v*) <L{u*,v*) + fi*^p{u*,v*) 

(4.3) = L{u*,v*,fi*) < snp{L{u,v,fi*) : {u,v) G C{U) x C{V)} 

= r, 

where the first inequality follows from the definition of Lagrange multiplier (^u* > 0) and the 
feasibility of {u*,v*) (i.e. tp{u*,v*) > 0). Using again the definition of Lagrange multiplier, 
we have /* = sup(„„)gj!f L(ti,t;,^*), so that equality holds throughout (j4.3p . This implies 
the equalities (jilD -dOD. 

Conversely, if {u*,v*) is feasible and (I4.ip - ()4.2p hold, we have from the definition of 
Lagrange multiplier, 

L{u*,v*) = L{u*,v*) + fi*i:{u*,v*) 

= L{u* ,v* , fi*) = max L{u,v, fi*) = L* , 

{u,v)ex 

so {u*,v*) is a global maximum. D 

Recall the definitions of the dual functional J in (I1.16P and the dual problem in (|1.17p . 
Note that J(^) may be equal to +oo for some /.i. In this case, we define the domain of J to 
be the set of fi for which J(/i) is finite: 

D = {fieR : J{fi) < +oo}. 

Regardless of the functional / and the constraint <j) of ths primal problem, the dual 
problem (I1.17P has a nice convexity property, as shown in the following lemma. 

Lemma 4.2. The domain D of the dual functional J is convex and J is convex over D. 

Proof. For any u,v,ii,p., and a € [0, 1], we have 

L{u, V, afi + (1 — a)/u) = aL(u, v, fi) + {1 — a)L{u, v, /2). 

Taking the supremum over all {u, v) E X, we obtain 

sup L{u, V, afi + (1 — a)/i) < a. sup L{u, v, fi) + {1 — a) sup L{u, v, p,), 

or equivalently 

J {a/j, + (1 — a)il) < aJ{ij) + (1 — a)J{jl). 
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Therefore if fi and ft belong to D, the same is true for a/i + (1 — a)p,, so D is convex. 
Furthermore, J is convex over D. D 

Another important property is that the optimal dual value 

J* = inf J(fi) 

is always an upper bound of the optimal primal value, as shown in the next lemma. 

Lemma 4.3. We have 

r < J*. 

Proof. For all /i > 0, and {u,v) E X with 'iIj{u,v) > 0, we have 

J{n) = sup L{u, V, /.f) 

iu,v)ex 

> I{u, v) + ^il^{u, v) > I{u, v), 



and therefore, 



J* = inf J{n) > sup I{u,v) = I*. 

/^>0 {u,v)£K 



D 



In the language of nonlinear programming [3l [10], if J* = I* we say that there is no 
duality gap; if J* > /* there is duality gap. Note that if there exists a Lagrange multiplier 
/i*, the above lemma (J* > /*) and the definition of Lagrange multiplier (/* = J{fJ-*) > J*) 
imply that there is no duality gap. 

The following is a sufficient condition for the existence of Lagrange multiplier, which is 
also a proof of Theorem 11.21 

Lemma 4.4. Let the assumptions ()1.18p - (|1.20p hold for the primal problem (J1.14p . Then 
there is no duality gap and there exists at least one Lagrange multiplier. 

Proof. Consider the subset of M^ given by 

A = {{z,w) : 3 {u, v) £ X such that 

ijj{u,v) > z, L{u,v) > w}. 

We first show that A is convex. Let {z,w) G A and {zjw) £ A he two different elements, 
we show that their convex combinations belong to A. 

The definition of A implies that for some {u,v) G X and {u,v) G X, we have 

L{u,v) > w, iIj{u,v) > z, 
I{u,v) > w, Tp{u,v) > z. 
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For any a € [0, 1], by the concavity of F in (jl.lSp . we obtain 

/ {au + (1 — a)u, av + (1 — a)v) > al{u, v) + (1 — a)I{u, v) 

> aw + (1 — a)w. 

By the convexity of ip in (|1.19p and noting that mfQ{f + h) > inf q / + inf q /i for any 
f,h£ C^{^), we obtain 

tp [au + (1 — a)u, av + {1 — a)v) > aip{u, f ) + (1 — a)'0('u, v) 

> az + {1 — a)z. 

Since the combination (au + (1 — a)u, av + (1 — a)v) £ X, the above equations imply that 
the convex combination of {z,uj) and {z,w), i.e. 

{az + (1 — a)z, aw + (1 — a)w) , 

belongs to A. This proves the convexity of A. 

We next observe that (0, /*) is not an interior point of A; otherwise, the point (0, /* + e) 
would belong to A for some small e > 0, contradicting the definition of /* as the optimal 
primal value. 

Therefore, there exists a supporting hyperplane passing through (0,/*) and containing 
A in one side. In particular, there exists a vector (/i, /3) ^ (0, 0) such that 

(4.4) f3r>nz + f3w, y{z,w)eA. 

We observe that if {z, w) G A, then (z,w — ^) £ A and {z — j,w) £ A for all 7 > 0. The 
inequality (|4.4p thus implies that 

(4.5) n>0, /3 > 0. 
We now claim that /? > 0. If not, /3 = and from ([0|) 

> HZ, y{z,w) G A. 

By the assumption (|1.20p . there exists a pair {u,v) G X such that 

ip{u,v) > 0. 

Since {'ip{u,v),I{u,v)) G A, we have 

> fitp{u,v), 

which in view of /x > in ()4.5p implies that fi = 0. This means, however, that (/i, /3) = (0, 0) 
arriving at a contradiction. Thus, we must have /3 > and by dividing if necessary the 
vector {fi, (3) by /3, we may assume that /3 = 1. Note that 

{ip{u,v),I{u,v)) £ A, y{u,v) £ X. 

Equation (14. 4p implies that 

(4.6) r >I{u,v)+iiij{u,v), y{u,v)£X. 
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Taking the supremum over {u, v) £ X and using the fact /i > 0, we obtain 

I* > sup L{u,v, fi) 

{u,v)ex 

= J{l^) > J*, 
where J* is the optimal dual value. By Lemma 14.31 we have the equalities hold above, 
namely /i is a Lagrange multiplier and there is no duality gap. D 

5. Examples and applications 

In this section we present some interesting examples and applications of the nonlinear 
optimization (fTTIl - ifO]) . 

5.1. Optimal transportation. Let U, V be two bounded domains in M", and c G C^{U x 
V) be a cost function. Let f,g be two positive densities supported on U,V, respectively, 
satisfying the mass balance condition 



(5.1) f= h- 

Ju Jv 
The optimal transportation problem is to find a measure preserving mapping Tq : U ^ V 

minimizing the cost functional 



C{T)= / c{x,T{x))f{x)dx 
Ju 

among all measure preserving mappings T from U to V. A mapping T :[/—>• 1/ is called 
measure preserving if 



h{T{x))f{x)dx = / h{y)g{y)dy 
lu Jv 

for any h G C{V). Denote the set of measure preserving mappings by T. 

Kantorovich introduced a dual functional 



(5.2) I{u,v)= / u{x)f{x)dx+ / v{y)g{y)dy 

Ju Jv 

over the set 

(5.3) K = {{u,v)eC{U)yiC{V) : u{x) + v{y)<c{x,y), 'ixeU,yeV]. 
Under suitable conditions, one can prove that 

inf C(T) = sup I{u,v). 

^^"^ iu,v)eK 

The reader is referred to [H El [TJ [HI \T9\ [22| [23] for further discussion on the optimal trans- 
portation problem. 

Note that ([52D-([53D is a hnear case of ([LT]) -([L2]1: 

F{x,y,u,v) = u{x)f{x) +v{y)g{y), d'y = dx®dy, 
(j){x,y,u,v) = u{x) +v{y) -c{x,y), 
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or equivalently, if{x, y,v) = v — c{x, y) in (|1.4p and (p*{x, y,u) = u — c{x, y) in (j3.2p . All the 
hypotheses in Theorem 1 1.1 1 are satisfied when the cost function c(x, y) satisfies the following 
conditions [T9]: 



(Al) For any x,p € M", there is a unique y £ M" such that Dxc{x,y) = p; and for any 

y,q £ M", there is a unique 2; G M" such that Dyc{x, y) = q. 
(A2) For any x, y G M", det D^^c / 0. 

The hypotheses (jl.3p - (jl.5p follow from the constructions of F, (j) together with / > 0, r^ > 0. 
The mass balance condition (|5.ip implies (jl.Sp . The condition (Al) implies both (HI) and 
(HI*), and (A2) implies (H2). 

Therefore, by Theorem II .11 we have the existence of potentials (u, v) and optimal mapping 

T. The existence of potentials in optimal transportation was previously proved in [HEllH]. 

By directly applying the formula (|2.16p . one obtains the optimal transportation equation 

f 
(5.4) I det [D'^u - d1^c]\ = |detc^,y|-. 

We remark that in the linear case (J5.2p - (j5.3p . both F in (jl.ip and (j) in (|l-2p are linear 
in t, s variables, that is a border situation of F being concave and being convex in t, s, 
simultaneously. 

There are numerous applicatons of the optimal transportation. Here we mention two 
important ones in geometric optics. In |25j . Xu-Jia Wang showed that the far field reflector 
problem is an optimal transportation problem, and so is a linear optimization problem. The 
associated cost function c{x, y) = — log(l — x • y), where x, y are points on the unit sphere 
§^. Later on in [12| Gutierrez and Huang showed that the far field refractor problem is also 
an optimal transporation problem. Let k be the refractor index from the initial media to 
the target media. Then the associated cost function c{x, y) = — log(l — kx ■ y) when k < 1; 
and c{x, y) = log(Ka; • y — 1) when k > 1, where x, y are points on the unit sphere S". In the 
following subsections, we will consider more general reflector and refractor problems and 
show that they are in the class of the nonlinear optimization (jl.ip - (jl.2p . 



5.2. Near field refiector problem with point source. In [16] we showed that the near 
field reflector problem is a nonlienar optimization. For the convenience of the reader, we 
summaries the arguments as follows. Assume that the light emits from the origin O and 
passes through C S*^ with a positive density / G L^(0). After being reflected from a 
surface F, the light will illuminate the target surface $7* in M*^"*"^ with a prescribed positive 
density g G L}{W). Assume the energy conservation condition 

(5.5) / / = / 9. 
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Represent the reflector T in polar coordinate system as 

Tp = {Xp{X) : X e n)}, 

where p is a positive function. Recall that [15], Tp is admissible if at each point Xp{X) G T 
there exists a supporting ellipsoid. Therefore, the radial function p satisfies 

(5.6) p(X) = inf ^^ ^^, X en, 

Yen'l-eipiY)){X,^)' 

where p is the focal function on Q* and e{p) = \J\ +p^/|yp — p/|y| is the eccentricity. 
Because there is an ellipsoid Ey^n^x supporting to V p for each y G fi*, we also have 



(5.7) p{Y) = sup p{X) 

xen 



l-eipiY)){X,^) 



Y £ n*. 



Note that in ()5.6p for each X £ Q the infimum is achieved at some Y £ Q* and in ()5.7p for 
each Y £ n* the supremum is achieved at some X £ n. 

The relations ()5.6p - ()5.7p between the radial and focal functions of a reflector T are anal- 
ogous to the classical relations between the radial and support functions for convex bodies, 
for example, see [21]. Inspired by that and [25], we set rj = 1/p. Then the pair (p, rf) satisfies 
the dual relation 

(5.8) p{X) = inf -. 

yen* „(v\ ( 1 _ ^(^ivwiv Y 



*r?(y)(l-e(7?(y))(X,^) 
r]{Y) = inf j -^ 

^e%(x)(i-6(^(y))(x,^)^ 

Similarly to |25j , we can now formulate the refiector problem to a nonlinear optimization 
(jl.ip ~ (jl.2p as follows. Let u = logp and v = log?/. Set the functional 



I{u,v)= [ F{X,Y,u,v) 
Jnxn* 



^^■^^ r ( (XY) \ 

= / f{X)u + g{Y) V + log(l - \l> =) , 

and the constraint set 

K = {{u,v)£ C(0) X C7(fi*) : 0(X, y, u, v) < 0} , 
with the constraint function 

(5.10) 4>{X,Y,u,v)=u^v^\oA\ ^^'J] |. 

We assume a further condition on domains il and n*: il* is contained in the cone 
CY = {tX : t>0,X £V} for a domain F C S" and 

(5.11) nnF = 0, 
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where ft and V denote the closure of il. and V. It imphes that there exists a small constant 
60 > 0, such that for any X £ Q,Y £ Q*, 

(5.12) - 1 < (X, ^) < 1 - 60. 

Under the assumption (j5.12p . one can verify that (jl.3p - (jl.5p are satisfied by (15. 9p and 
(jS.lOp . The energy conservation condition (15. 5p implies (ll.Sp . To show the condition (HI), 
let {u,v) G if be a dual maximizing pair of /, (see Lemma l2.ip . If there holds 

MX,Y,v{Y)) = -Du{X), ip{X,Y,v{Y)) = -u{X), 

at X £ il.,Y £ il.* , it was proved |16j that Y is the target point of light emitting along X, 

reflected at Xe"' ' £ T with unit normal 

_ {Du,0) -il + Du-x)X 

^~ y/l + \Du\^-{Du-xy' 

By the reflection law, Y,. = X — 2{X,j)-y, 

Y = Xe'^W +Yr\Y - Xe"(^)| 

is uniquely determined. It implies that {ipx, ^)i^: "; ") is one-to-one in i}* x v{i}*) for each 
X £Q, see Remark [221 

As a consequent of Theorem II. H we have 



Corollary 5.1. Assume that f,g satisfy (j5.5p . Suppose thatQ andQ* satisfy (|5.1ip . Then 
there is a dual maximizing pair (u, v) £ K satisfing 

I{u,v) = sup I{u,v), 

{u,v)£K 
where I{u,v) is in (j5.9p . and the constraint <j) is in (j5.10p . Moreover, p = e^ is a solution of 
the reflector problem with given densities (^i, /) and {Q*,g). (Note that the solutions need 
to be understood as generalized solutions.) 

By directly applying the formula ()2.16p . we also obtain the PDE in the near field case 
|16| . which was previously obtained by Karakhanyan and Wang in [15]. Assume that O* is 
given implicitly by 

(5.13) n* = {Z £ IR"+^ : ^{Z) = 0}. 

Suppose that fl is a subset of upper unit sphere §!f: = S" n {x„+i > 0}. Let X = (x, x„+i) 
be a parameterization of fi, where Xn+i = y/l — \x\^ =: uj{x), and x = [xi, ■ ■ ■ ,x„). For 
simplification, we define some auxiliary functions 

(5.14) a = \Dp\^-{p + Dp-xf, 

(5.15) b = \Dp\^ + p^-{Dp-xf, 

(5.16) t = ^""+^~^"+\ /3 



pxn+i ' "^ {Y-Xp)-V'^' 
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and denote the matrix 



(5.17) N = {Mj}, N^J = 6ij + ^_^ 

By computing in the local orthonormal frame, we obtain the equation as follows 

Corollary 5.2. The function p is a solution of 



(5.18) 



det 



n 2 a(l- 1) ' 

D^p - -Dp ®Dp- \^ ' M 
p 2tp 



^n+l 



t^bp 



f 



2"/)2n.+1^2^|V^|' 



Note that the matrix in equation ()5.18p has a different sign to that in p!5], since we 
calculate the absolute value of the determinant. 

5.3. Near field reflector problem with parallel source. The ideal reflection system can 
be described as follows: a parallel light emits from C M" x {0} along e„+i = (0, • • • ,0,1) 
with a positive density / G L^{Q). After being reflected by the surface T in M'""^"^, it 
will illuminate the target domain Q,* C M" x {0} with the prescribed density g G L^(il*). 
Assume the energy conservation condition 



(5.19) / / = / g. 

Jn Jn* 

We represent the reflector T as graph u\^, namely 

Tu = {{x,u{x)) : X gQ}, 

where n is a positive function. 

Consider the "inverse" paraboloid with focus at y £ Q* and the axial direction —Cn+i- 
The reflection property of these paraboloids is that: all the incident light from parallel 
source along the direction e„+i will be reflected by the inverse paraboloids to the focus 
points y. Such an inverse paraboloid can be represented by Tp = {{x,p{x)) : x G M"} with 

(5.20) p(x) = py{x) = — - -\x - yl"^ , 

where f > is a constant depending on y such that p{x) > in il. 

If we regard v = v{y) as a function on Q* , we then have a family of paraboloids {Tpy : y £ 
0,*}. It is natural to construct the reflector r„ by the envelope of the family of paraboloids 
{Tp, : y G n*}. 

In an ideal system, at each point (x,u{x)) G r^ on an admissible reflector Tu there exists 
a supporting paraboloid, namely for some y £ W 

^^"^^ = 24)"^'"^"^''' """"^ n(^')<^-^k'-yp Vx'Ga 
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Hence, the defining function u and the dual function v satisfy the dual relation: 

(5.21) u{x) = mi\-^-^\x-yA, x e Q, 

y^n* [2v[y) 2 J 



(5.22) .(,)= inf Z!i±Y?+fH^ I yen*. 

xen \x — y|^ I 

Note that in (j5.2ip . if for x e il. the infimum is achieved at some y = T{x) € il.* , then in 
(|5.22p at y = T{x) the infimum is achieved at x € il. By ()5.2ip one has 

Du{x) = —v{x — y), 

and thus 

\Du{x)\'^ = l-2vu. 
Since n > 0,t; > 0, we have the natural restriction \Du\ < 1 on ri. 

Based on ()5.2ip - ()5.22p we can now formulate this problem to a nonlinear optimization 
problem ()l.ip - ()1.2p . Set the functional 



I{u,v) = j F{x,y,u,v) 
(5.23) * 



f{x)u{x) + g{y) \-\x-y 



and the constraint set 



V, '2 _ ]_ 

2v 



/C = {{u,v) E C{n) X C(0*) : (t){x,y,u,v) < 0}, 

with the constraint function 

1 V 
(5.24) (t>{x,y,u,v) =u- — + -\x -y\'^. 

2v 2 

Note that can be written as (|1.4p with (p{x,y,v) = —i^ + ||x — yp. By differentiating, 

we have 



^s = Trpi ^ '^\^ ~ y\ ' ^xs = x-y, etc. 



(5.25) ifj; = {x-y)v, (fy = -{x-y)v, i.pxx = vl, (pr,y = -vl, 

1 1 

2"^ ^2' 

From Lemma |2. 1 1 there exists a dual maximizing pair (n, v) £ K oi I. From (|5.25p one can 
see that {(px,^){x, •, •) is one-to-one in Q* x i;(il*) for each x E il, Remark 12. 2i The energy 
conservation condition (|5.19p implies (jl.Sp . Alternatively, one can directly verify that the 
optimal mapping T associated to {u,v), determined by (|2.6p . is equal to the reflection 
mapping. This implies the condition (HI) due to the reflection law. 

Proposition 5.1. Let {u,v) £ K be a dual maximizing pair of I, \Du\ < 1. The associated 
optimal m,apping T determined by (|2.6p is equal to the reflection mapping Tr obtained by 
the reflection law. 
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Proof. From (j2.6p and (j5.25p . at {x,y) = {x,T{x)) we have 

Du 

V 



(5.26) 



x-y 



From ([LTD, u - ^ + ^\x - y\^ = 0, thus 



(5.27) 



1- \Du\ 
2u 



Note that u > 0,v > and \Du\ < 1. Combining (15.26J) and (|5.27p . we obtain 
(5.28) '"^" 



x-y 



l-\Du\^' 
This implies the optimal mapping T is given by 

2uDu 



T{x) = X + 



1- \Du 



2' 



Next, we calculate the reflection mapping T^. At point [x,u{x)) G F^, from direct 
calculations the unit normal is 

^l + \Du\^' 
By the reflection law, we have the reflected direction 

{2Du,\Du\'^ - 1) 



(5.29) 



Yr = en+i - 2(e„+i,7)7 



\Du\^ + l 



On the other hand, since the reflected ray meets the hyperplane R" x {0}, we have 

{y-x,-u) 



(5.30) Yr 

From (fOOjl - ljOO]) we obtain that 
(5.31) 



^\y — xp + u'^ 



r^ , N 2uDu „, , 

y = Tr{x) =X + — - = T X . 

1 — \Uur 



U 



Therefore, as a consequence of Theorem II. H we have 



Corollary 5.3. Assume that f,g satisfy the energy conservation condition (|5.19p . Then 
there is a dual maximizing pair {u,v) G K of I, where I,K are in ()5.23p . ()5.24p . Moreover, 
u is a (generalized) solution of the reflector problem with given densities (^i, /) and {Q*,g). 

We now derive the equation for the potential function u by using the formula (j2.16p . 



|det [D'^u + vl] 



(5.32) 



det 



1 

fs 



(x -y)0{x- y) 



1 \x — y\ 



f 
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where in the second step we used the formula det [I + £, <i$ ij] = 1 + ^ ■ rj for any vectors 
^,r] € W. Combining ([OSjl - fjOH]) . we obtain that 

l|a;-yP 1, ,9 l + l-Dup, „ 
(5-33) c^s = 2^0^ + 2^ - 2/1 = ^^^1^ " 2/1 • 

Therefore, the equation (|5.32p becomes 



(5.34) det -^^ ' ^ '^-^1 ^ ^' l^""' ^ ^ 



D^u+ ^ ^/ 



2u 



(2u)"(l + |L»u|2)g' 



The equation (|5.34p was also obtained in [18] by directly computing the Jacobian of the 
reflection mapping. Denote the matrix 



|2 



A{u,Du)= ^ j^^' I. 

It satisfies the (A3) condition in [19] without the orthogonal restriction, provided that 
u > 0. The regularity of (j5.34p follows from [18^ I19j . Indeed the considerations in [19] 
stemmed from the treatment of the reflector antenna problem by Wang in [2l] , which can 
be represented as an optimal transportation problem on the sphere §" with the cost function 
c{x,y) = -log(l -x-y). 

5.4. Near field refractor problem with point source. The near field refractor problem 
has been studied by Gutierrez and Huang [13] , Suppose the light emits from the origin 
surrounded by medium I with positive intensity f{X) for X £ Q, where 0, C S". There is 
a surface TZ, separates two homogeneous and isotropic media I and II, such that all rays 
refracted by 7^ into medium II illuminate a target hypersurface ri* in R'"+^ with positive 
intensity g on Q* . Assume that /, g satisfy the energy conservation condition 



(5.35) / / = / 9- 

Let ni,n2 be the indices of refraction of media I, II, respectively, and k = n2/ni. When 
K < 1, the refracted rays tend to bent away from the normal, when k > 1, the refracted 
rays tend to bent towards the normal. 

There is a special interface surface S between media I and II, called Cartesian oval [13j . 
that refracts all rays emitting from the origin O into the point Y. In the polar coordinates, 
represent 

S = {XpoiX) : X G S"}. 

In the case k < 1, by the Snell law of refraction one has 



(5,36) ..m-Mxr.rt^ '^-^-^-f'-.v^^^^^ . 

where p is the focal parameter, k\Y\ < p < \Y\ and A(t) = {p — K^t)'^ — (1 — K^)(p^ — K^|yp) 
for t G M. For non-degenerate Cartesian ovals, there are some physical constraints for 
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refraction 

(5.37) k\Y\<p<\Y\, Po<P, andp<X-Y, VX £ n. 

We refer the readers to [13] for more physical interpretations and detailed calculations. 

If we regard p = p{Y) as a focal function on il*, we then have a family of Cartesian ovals. 
Represent the surface TZ in polar coordinate system as 

Up = {Xp{X) : X €n}, 

where p is a positive function. Recall that |13j . 7^ is a near field refractor if at each point 
Xp{X) € TZ there exists a supporting Cartesian oval, i.e., for some Y £ il.* , p{X') < 
Po{X' ,Y,p{Y)) for all X' £ Q. with equality holds at X' = X. Therefore, the radial function 
p satisfies 



(5.38) piX) = mf ^--'^■y)-/WrYl ^ ^ ^ 

yen* l-K? 

From the energy conservation, for each y G il* there is an oval pQ{-,Y,p{Y)) supporting to 
TZp. We also have 

(5.39) p{Y) = sup (1 - K^)p{X) + k'^X ■ Y + ^J /\{X ■ Y), Yen*. 

xen 

The above relations are analogous to (|5.6p - (|5.7p . By setting t] = 1/p, the pair (p, rj) satisfies 
the dual relation 

(5.40) .(X) = i„f '-"-"^■f^f^^^ ll. 

Y&Q* r]{l - K^) 



Similarly to (j5.9p we now formulate the refractor problem to the following nonlinear op- 
timization, which is more complicated than (j5.9p . Let u = logp and v = \ogrj. Set the 
functional 

Iiu,v)= [ F{X,Y,u,v) 

and the constraint set 

K = {{u, v) G C{n) X C{n*) : (p{X, Y, u, v) < 0}, 
with the constraint function 

(5.43) (/>(X, Y,u,v) =u + v + log 



1 - e«(K2X • Y + ^A{X-Y))J 

As in [13] , we make the following assumptions on and 0* , which are due to the physical 
constraints for refraction (15.371): 
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(Rl) There exists r with < r < 1 - «; such that X -Y > {k + t)\Y\ for ah x G H and 

ally eQ*. 

(R2) Let < ro < Y^dist(0,n*) and consider the cone in M"+^ 

Qro={tX : X en,0 <t< ro}. 

For each ^ G M" and for each X G Q^o we assume that 0, D {X + t^ : t > 0} 
contains at most one point. That is, for each X G Qr^ each ray emanating from X 
intersects O at most in one point. 

Note that from (f06]) 

(5.44) A(X • Y) = kV - 2k2(^ • Y)p + k^{X ■ Yf + ^^(1 - k^)\Y\'^. 

Since p < X -Y hy (j5.37p . A(X ■ Y) is decreasing in p. Since p = e~'", A{X ■ Y) is increasing 
in V. Hence (|5.42p - (|5.43p satisfy the assumption (|1.3p and (|1.5p . where the constant in (|1.5p 
depends on K,Q,fl* and the assumptions (Rl) and (R2). See Remark 4.2 in [13] for more 
physical interpretations of (Rl) and (R2). 

By the proof of Theorem 11.11 and Remark 12. II we have the following existence result in the 
near field refractor problem, which was previously obtained in [13] . Note that the solutions 
need to be understood as generalized solutions. 



Corollary 5.4. Assume that f,g satisfy (j5.35p and (R1)-(R2). Then given Pq G Qtq 
with < I Pol ^ ( TTk ) '"Oi there exists a dual maximizing pair {u,v) G K such that 
u{Po/\Po\) = log |Po|- Moreover, p = e^ is a solution of the refractor problem such that IZp 
passes through the point Pq. 

Remark 5.1. In the case k > 1, by the physical constraints for \Y\ < p < k\Y\ the refracting 
piece of Cartesian oval is given by 

with 



(5.45) p.(xr.p) = "-'^' ^ ^ -"' - ^^<"^^ ■ \-'^'' - '-' - i)(-VI^ -J-^). 

We have similar existence results by replacing the assumptions (Rl), (R2) by the following 
ones: 



(R3) mij^^^y^^* X • ^ > i + T, for some < r < 1 - i. 



.2^2 



(R4) Let < ro < 4^(l-IiV2 infyen* 1^1 and consider the cone i 



m 



pn+l 



Qro = {tX : X £n,0<t< ro}. 
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For each ^ G M" and for each X G Q^o we assume that $7 n {X + t^ : t > 0} 
contains at most one point. That is, for each X G Qro each ray emanating from X 
intersects 0, at most in one point. 

5.5. Near field refractor problem with parallel source. Recently, Gutierrez and 
Tournier studied the parallel refractor problem [1^, which can be described as follows. 
Suppose that a parallel light emits from Q C ffi"" x {0} along e„+i = (0, ••• ,0, 1) with 
positive intensity / G L^(0), Q* is a hypersurface in M""'"^, which is referred to as the target 
domain. Suppose that fJ and O* are surrounded by two homogeneous and isotropic media 
/ and //, respectively. One seeks an optical surface TZ interface between media / and //, 
such that all rays refracted by 7?. into medium // are received at the surface Q* , and the 
prescribed radiation intensity received at each point 1" G O* is g{Y). Assume the energy 
conservation condition 



(5.46) / / = / 9- 



Let ni,n2 be the indices of refraction of media I, II, respectively, and k = ni/n2. We 
assume that media // is denser than media /, that is, «: < 1. The case when k > 1 can 
be treated in a similar way but the geometry of surface changes |121 ^^ . For simplicity, we 
assume that i7* C {Vn+i = ^} for a constant /i > 0, and denote Y = (y, h) for points on VL* 
and X = (x, 0) for points on 0. 

Consider the lower part of "inverse" ellipsoid of revolution with focus at y G il* and 
the axial direction — e^+i- It has the uniform refracting property, namely all rays from the 
parallel source il along e„+i will be refracted to the focus points y. Explicitly, it is the 
graph of the function 



Kv / f 2 \x — yp 



where w is a constant satisfying ^ ~"^ ' < v < ^ ~'^ 3 — ^. The function py^^ is defined on 
the baU B,^i^^^^{y). 

If we regard v = v{y) as a function on il.* , we then have a family of "inverse" ellipsoids. 
Represent the refractor F as graph u\n for n > 0, namely 

Fu = {{x,u{x)) : X G 0}. 
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In an ideal system, at each point {x,u{x)) G r„ on an admissible refractor r„ there exists 
a supporting ellipsoid, i.e., for some y £ ^* 






Kv{y) / v{y)'^ \x' — yP 






Similarly we can formulate this problem to a nonlinear optimization problem ()l.ip - ()1.2p . 
Set the functional 



I{u,v) = I F{x,y,u,v) 
Jnxn* 
(5.48) 



Jnxn* \ 1 - K^ 



i^v{y) / v{yy \x - y\ 



2 



(1 - k2)2 1 - k2 

and the constraint set 

K = {{u, v) G C{n) X C{n*) : <p{X, Y, u, v) < 0}, 
with the constraint function 



(5.49) HX,Y,u,v)=u + ^ + J^^-\-^-h. 



As in [14] we need the following assumptions on the relative position of Q, and Q,* , which 
are due to the physical constraints for refraction: 

(A) There exists < S < 1 such that C -B^/iVl - k? / n{y) for ah y G il*. 

(B) Set M = /i((l + nf/n^ - 1). Assume that for ah x G x [-M, 0] and for ah 7 G S", 
the ray {x + ^7 : i > 0} intersects Vl* in at most one point. 

The first condition is equivalent to the assumption that there exists < /3 < 1 such that 
(-e„+i, ^5^) > /3 for ah y G Jl* and 2; G Jl il4j. 

As in the previous example one can verify the hypotheses of Theorem ll.il Therefore, by 
the proof of Theorem 11.11 and Remark 12. II we have the existence result: 

Corollary 5.5. Assume that f,g satisfy ()5.46p and (A)~(B). Then for xq G and t < —fi 
there exists a parallel refractor u satisfying u{xq) = t. 

Note that the solutions need to be understood as generalized solutions as before. This 
existence result was previously obtained in [H], where they first consider the discrete case 
when the target is a set of points, then use an approximation to obtain the existence in the 
general case. 
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Remark 5.2. In addition to the examples arising in reflectors and refractors, there are many 
other nonlinear optimization problems with potentials. For example, one can perturb the 
linear optimization problem, such as the optimal transportation, to get a nonlinear one. 
Moreover, similarly to [2^ one can show that the objective functional of any solvable linear 
optimization problem can be perturbed by a differentiable, convex or Lipschitz continuous 
nonlinear functional in such a way that (i) a solution of the original linear problem is a 
local or global solution of the perturbed nonlinear problem; (ii) each global solution of the 
perturbed nonlinear problem is also a solution of the linear problem. 
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